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ABSTRACT 

Context. A Band function has become the standard spectral function used to describe the prompt emission spectra of gamma-ray 
bursts (GRBs). However, deviations from this function have previously been observed in GRBs detected by BATSE and in individual 
GRBs from the Fermi era. 

Aims. We present a systematic and rigorous search for spectral deviations from a Band function at low energies in a sample of the 
first two years of high fluence, long bursts detected by the Fermi Gamma-Ray Burst Monitor (GBM). The sample contains 45 bursts 
with a fluence greater than 2xl0~ 5 erg / cm 2 (10 - 1000 keV). 

Methods. An extrapolated fit method is used to search for low-energy spectral anomalies, whereby a Band function is fit above 
a variable low-energy threshold and then the best fit function is extrapolated to lower energy data. Deviations are quantified by 
examining residuals derived from the extrapolated function and the data and their significance is determined via comprehensive 
simulations which account for the instrument response. This method was employed for both time-integrated burst spectra and time- 
resolved bins defined by a signal to noise ratio of 25 a and 50 a. 

Results. Significant deviations are evident in 3 bursts (GRB 081215A, GRB 090424 and GRB 090902B) in the time -integrated sample 
(~ 7%) and 5 bursts (GRB 090323, GRB 090424, GRB 090820, GRB 090902B and GRB 090926A) in the time-resolved sample (~ 
11%). 

Conclusions. The advantage of the systematic, blind search analysis is that it can demonstrate the requirement for an additional 
spectral component without any prior knowledge of the nature of that extra component. Deviations are found in a large fraction of 
high fluence GRBs; fainter GRBs may not have sufficient statistics for deviations to be found using this method. 

Key words. Gamma-ray burst: general - Methods: data analysis - Techniques: spectroscopic 



1. Introduction 

Gamma-ray bursts (GRBs) are the most luminous events in the 
universe and can briefly be summarised as having high-energy 
prompt emission followed by a multi-wavelength fading after- 
glow (e.g., |Vedrenne & Atteia||2009) |Kann et al.|[20Tl) . The 
isotro pic energy produc ed by a typical GRB is E iso ~ 10 51 erg 
(e.g., Frail et al.pOOl and in some cases reaches up to 10 



erg (e .g., |Greiner et al.|2009||McBreen et al.|2010||Cenko et al.| 
201 1) . The prompt emission of GRBs has been detected over 



a wide spectral range from keV to GeV energies and is gener- 
ally well modelled by one or a combination of the following: a 



smoothly broken power-law (e.g. Band et 


al.|1993 Abdoetal. 


2009b), a quasi-thermal component (e.g., 


Preece 2000 Guiriec 


et al.|2011||Ryde et al.|201 l||Axelsson et a 


.|2012|l, an extra non- 


thermal power-law component extending 


to high energies (e.g. 



|Gonzalez et al.||2003| |Kaneko et al.| 


2008 Abdo et al. 2009a) 


or a cut-off in the MeV regime (e.g., 


Ackermann et al.|2012b). 



The lightcurves of GRBs are highly variable and a number of 
studies into their temporal properties have been performed (e.g., 
IQuilligan et al.|2002[|Hakkila & Preece|201 l||Bhat et al.|2012) . 

Our work assumes that a Band function ( |Band et al.|1993| l is 
the best fit function for GRB spectra in the GBM energy range 
and attempts to quantify the number of bursts that deviate from 
this function. In general, a Band function can be constrained bet- 
ter at lower energies compared to higher energies as more counts 
are observed at lower energies. This provides strict limits on the 
a parameter and makes it possible to search for deviations to the 
fit function. The Band model is defined as 



N(E) 



if E < E r 



U(^r«po-«)(Ay if e> Ec 



where 



E c = (a-0) 



^ peak 

2 + a 
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and A is the amplitude in photons s _1 cirT 2 keV -1 , 
a is the low-energy power-law index, E pe ak is the v F v 
peak energy in keV and (3 is the high-energy power-law index. 

The emission mechanisms for GRBs have generally been de- 
scribed as thermal or non-thermal with the thermal emission 
characterised by the radiation emitted from a cooling plasma 
(e.g. Meszaros & Rees 2000). However, thermal emission need 
not take the form of a standard Planckian model but may take a 
more complex form depending on several factors including the 
relative leptonic and hadronic populations, viewing angles and 
source region (e.g., Ryde et al. 201 Non-thermal emission can 



Table 1. Comparison of properties of the GBM Nals with the 
BATSE SDs. 



take the form of synchrotron emission from the extreme mag- 
netic fields needed to create a GRB (e.g. Meszaros 2002 ). Strong 
magnetic fields may also break and reconnect releasing energy 
in the form of gamma-ray photons (e.g. |Zhang & Yan 201 1) . 
Although a Band function is not a physical model, it successfully 
models a large number of GRBs and it is useful to investigate 
the number of events which require additional paramaters. If an 
additional blackbody is observed, the temperature of the photo- 
sphere can be deduced ( |Rees & Meszaros 2005 ) whereas an un- 
derlying power-law can provide constraints on emission mecha- 
nisms (e.g. |Asano et al.|2009||Ryde et al.|2010[) and the E xtra- 
galactic Background Light (EBL) ( |Ackermann et al.|201 1] >. An 
overview of the observational and physical interpretations of 
the results from the Fermi era are presented by Bhat & Guiriec 
pOTT) . 

A number of papers discuss interesting bursts with unusual 
spectral behaviour in the keV range (e .g., Guiriec et al. |2011 



Ryde et al.|2011| |Axelsson et al.||2012[ >. It is important to note 



that there is a bias in the literature whereby interesting bursts are 
published more than bursts that conform to the standard GRB 
models. Another bias that can skew the apparent number of in- 
teresting bursts is that only bursts that look atypical initially are 
investigated further. A rigorous investigation must be carried out 
on a sample of GRBs in a systematic way to obtain the true frac- 
tion of GRBs with observable additional features. Any investi- 
gation must also be performed blindly to decrease the risk of 
interval selection effects which could bias the study. Here we 
present a systematic search for deviations from a Band function 
in the highest fluence bursts in the first catalogs from the Fermi 
Gamma-ray Burst Monitor ( Paciesas et al.|2012 Goldstein et al. 
|20T2] l. 

I. 1. Fermi Gamma-ray Burst Monitor 

The Fermi Gamma-ray Space Telescope, launched on 2008 June 

II, has an energy range spanning several decades (~8 keV to 
~ 300 GeV) and is ideal to explore the low-energy regime 
of GRBs. Fermi consists of two instruments, the Large Area 
Telescope (LAT) operating between ~ 20 MeV to ~ 300 GeV 
( |Atwoodet al.|2 009 ) and the Gamm a-Ray Burst Monitor (GBM) 
operating between 8 keV - 40 MeV ( |Meegan et al.|2009) >. 

GBM consists of two types of detectors - twelve Sodium 
Iodide (Nal) scintillating crystals operating between 8 - 1000 
keV and two Bismuth Germanate (BGO) scintillating crystals 
operating between 0.15 - 40 MeV. The Nal detectors are ar- 
ranged in clusters of three around the edges of the satellite and 
the BGOs are located on opposing sides of the satellite aligned 
perpendicular to the LAT boresight. As the detectors have no ac- 
tive shield and are uncollimated, GBM observes the entire unoc- 
culted sky. Further details can be found in Meegan et al. (2009 ). 

Both Nal and BGO detectors collect counts in 4096 channels 
compressed into 128 energy channels, with boundaries spaced 





BATSE SD 1 


GBM Nal 2 


Material 


Nal 


Nal 


Number 


8 


12 


Area 


126 cm 2 


126 cm 2 


Thickness 


7.62 cm 


1.27 cm 


Energy Range 


30 keV - 10 MeV* 


8 keV - 1 MeV 


Spectral Binning 






8 - 20 keV 


1 bin 


~ 12 



Notes. (*) This energy range varied depending on the gain of the detec- 
tors where the lower end of the range could be as low as 5 keV for a 
high gain setting iPr eece et al.| 1 996 ). 

References. (l JPreece et ^T) ( [T996l >, (2) |Meegan eTaT] ( (2009^ 



pseudo-logarithmically across the energy ranges of each detec- 
tor. The lowest and highest energy channels in each detector 
are generally ignored in the analysis because of uncertainties 
in the instrument response. Extensive ground calibration was 
carried out on the detectors pre-launch (Bissaldi et al. 2009) 
and in flight by comparison to other instruments (INTEGRAL- 
ISGRI : |rTerney et al.| ( [20TT|>, INTEGRAL -SPI: |von Kienlin et al.| 
pD09) i and Swi^B AT ^tamatikos| ( [20091 ) . These calibration re- 
sults are consistent and do not show any major changes from 
the on-ground calibration. Additional calibration work has also 
been carried out using the Crab Nebula with Earth occultation 
techniques ( Wilson-Hodge et al.|2012| l, nuclear lines from solar 
flares (Ackermann et al. 2012a| l, positron/electron annihilation 
lines (Bri ggs et al.|201l" l and spectral analysis of Soft Gamma- 
ray Repeaters (SGRs) ( |Lin et al.|2012| i. 



1.2. Fermi GBM and BATSE/CGRO Comparison 

A comprehensive study examining deviations at the lower end of 
the spectrum has not been carried out since Preece et al. ( 1996) 
investigated time-integrated GRB spectra using the spectroscopy 
detectors (SDs) on the Burst And Transient Source Experiment 
(BATSE) on board the Compton Gamma Ray Observatory 
(CGRO) ( |Gehrels et al.|1992| ). The results of this study showed 
that ~ 14% of a sample of 86 GRBs contained significant low- 
energy excesses; no significant deficits were observed in the 
sample. A comparison between the effective areas of BATSE and 
GBM is given in Figure[T] The GBM Nal detectors are of similar 
specifications to the BATSE SDs (see Table [T}. As the BATSE 
SD detectors were significantly thicker than the GBM Nal de- 
tectors, the BATSE SDs could detect higher energy photons than 
the GBM Nals. This is compensated for on GBM by using the 
BGO detectors for high-energy constraints. 

The relative effective areas for each detector are presented 
in Figure [T] showing that in the crucial energy range of interest 
(below ~ 30 keV), the Nals have a greater effective area than the 
SD detectors. Additionally, while Preece et al. (1996) only used 



a single SD detector, multiple Nal detectors were used in our 
analysis in all but 3 GRBs (see §[2). GBM also has the advantage 
of greater effective area at higher energies using the BGOs. This 
provides better constraints on /3 and E pea k which help constrain 
the overall spectral model applied to the data. 

To probe the lower energies using the BATSE SDs, two data 
types were required. The first data type was the Spectroscopy 
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Fig. 1. Plot comparing the efFective areas of GBM v BATSE per 
detector. The total area exposed to a particular sky direction de- 
pends on the orientation of each detector and number of detec- 
tors exposed. At low energies the GBM Nal detectors have more 
effective area than the BATSE SDs and at high energies the GBM 
BGO detectors have more effective area that the BATSE LADs 
( Meegan et al.|2009 1. The dark lines represent the photopeak ef- 
fective area where the detected energy is the same (within the 
energy resolution) as the incident energy and the lighter lines 
represent the total effective area which includes the photopeak 
plus the cases where the instrument detects only part of the inci- 
dent photon energy. (A color version of this figure is available in 
the online journal.) 



or more Nals were used for 42 out of the 45 GRBs. The bright- 
est BGO detector was used in all cases to constrain the spectral 
model at high energies. 

As the signal to noise ratio (S/N) in these bursts is quite 
high, the influence of background fluctuations has a less statis- 
tically significant effect during the spectral fitting process com- 
pared to weaker bursts. Individual detectors were checked for 
blockages whereby the source photons passed through part of 
the spacecraft before entering a detector. This can cause photons 
to be scattered to different energies and lower energy photons 
to be absorbed beyond the ability of the detector response ma- 
trices (DRMs) to accurately reconstruct the photons incident on 
the entire spacecraft. An automated tool was initially used to 
check for line of sight blockages between the GRB and a detec- 
tor. Further manual analysis was performed to remove any addi- 
tional blocked detectors. Detectors that are blocked have a not- 
icable deficiency of low-energy counts which is not consistent 
with other detectors that have similar source angles. Detectors 
which were not blocked and had an acceptable source angle were 
defined as 'good' and used in the subsequent analysis. 

Figure [2] shows the spectral properties of this sub-sample 
(a and E pe! ^) re lative to the overall GBM spectral catalog 
( |Goldstein et al.|[2012[ ). The distributions show that the sub- 
sample selected here displays similar characteristics to the entire 
GBM spectral catalogue. 



3. Method 



Discriminator data channels (DISCSP) which contained four in- 
tegral channels over the entire energy range ~ 5 keV - 2 MeV 
(for the highest gain setting). The upper edge of the lowest data 
channel (DISCSP 1) was set to the Lower- Level Discriminator 
(LLD) threshold of the Spectroscopy High-Energy Resolution 
Burst (SHERB) data (~ 10 keV for the highest gain setting). The 
SHERB data comprised 256 spectral energy channels between 
the LLD and ~ 2 MeV. By using a joint fit between the SHERB 
data and the single DISCSP 1 data point, |Preece et al.| (1996) 
determined the deviation of this single data point to the model. 
GBM has the advantage of a single standard data type from low 
to high energies. This reduces the chance of any incongruity be- 
tween data sets. The data type also has better spectral resolution 
in the region of interest (see Table [TJ so trends in the data can be 
more easily distinguished. GBM also has the advantage of bet- 
ter source locations provided by other instruments such as Swift 
and the LAT which improves the ability to accurately model the 
instrument response. Out of the sample of GRBs analysed in this 
work, 47 % have a sub-degree localisation. 



2. Sample Selection 

The sample was drawn from the first 2 years of triggered GBM 
GRBs (U th July 2008 - W h July 2010). Bursts above a fluence 
of 2 x 10~ 5 erg cirT 2 (10 - 1000 keV) were selected so that 
the spectral parameters could be constrained in the fitting pro- 
cess (see however § 3.5 1. This gave a sample of 45 bursts, which 
formed the brightest 9% of th e 491 GRBs in the first GBM GRB 
catalog ( |Paciesaset al.|2012| >. In addition, 36 out of the 45 GRBs 
(80%) have significant emission in at least one BGO detector as 
defined by Bissaldi et al.|(|20lT). Only Nal detectors with source 



The sample was analysed using two methods. The first compares 
the data to the extrapolated fit of the function in the low-energy 
regime and the second compares the change in the spectral index 
a when the energy range of the data is shortened. These two 
methods are used in combination to identify deviations from a 
simple power-law at lower energies. 

For an additional component to be significantly detected in 
the spectrum of a burst, it must either be continuously present 
throughout the entire burst or very strongly present in certain 
sections of the burst. An additional component may be present 
in certain sections of a GRB but not be intense enough to sig- 
nificantly alter the overall spectrum of the burst. Therefore the 
spectral fitting was performed on the spectrum of the entire burst 
(time-integrated) and on significant time slices (time-resolved) 
of the data. All spectral fitting was performed using the RMFIT 
software package [|] and minor modifications were made to auto- 
mate the fitting process. 

3.1. Time-Integrated Single Fit 

CSPEC data were used for each 'good' detector over the full en- 
ergy range of GBM (8 keV - 40 MeV). These data have an energy 
resolution of 128 channels and a temporal resolution of 4 s pre- 
trigger changing to 1 s post-trigger until TO + 600 s. Using the 
lightcurve, a time region encompassing the main emission phase 
of the GRB was selected by eye. A polynomial background was 
then fit to each GRB lightcurve by selecting background regions 
before and after the prompt phase. Multiple response matrices 
were used in the fitting process to account for spacecraft slew- 
ing relative to the source. In this process a new response is made 
when the spacecraft slews by more than 2 degrees, which can be 
important in long bursts, and is especially important if the space- 



angles less than 60° were used in the spectral analysis to limit the 
effect of the uncertainties in the off-axis detector response. Two 



http://fermi.gsfc.nasa.gov/ssc/data/analysis/user/ 
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Fig. 2. Histograms comparing a) a and b) E pe ak in the high fluence sample from our work to a sample obtained from 



Goldstein et al. 



craft executes an autonomous repoint. A similar method was also 
employed by Goldstein et al. ( 2012| l. 

An initial time-integrated spectral fit was performed using a 
Band function for each GRB in the sample. Fit residuals were de- 
termined by the number of standard deviations that separate the 
data from a Band function. The residuals were summed between 
8 keV and a variable Low-Energy Threshold (LET) in order to 
search for excesses or deficits. The LET was set at 15, 20, 25, 
30, 50 and 100 keV. However, this technique is not optimal as 
the fitting algorithm forces the function through any deviations 
that are present and minimises any deviations in any particular 
part of the spectrum. To reduce this issue, the LET was used as 
a lower energy bound for spectral fitting and the data points be- 
low the LET were compared to an extrapolated version of the 
function. This extrapolated fit is defined in the next section. 

Although the sample was defined to be composed of high 
fluence GRBs, this does not necessarily guarantee a large num- 
ber of counts in all energy channels. Due to the potential low 
count rate / bin ratio, the fit (minimization) was performed using 
Castor statistics. This is similar to Cash-statistics (Cash 1979) . 
Castor statistics assume Poisson uncertainties per bin compared 
to x 2 statistics which assumes Gaussian uncertainties per bin. 
These are equivalent in the high count regime but diverge in the 
low count regime (less than ~ 10 counts / bin). A similar min- 



imisation method was used by Goldstein et al. ( 2012] >. 



3.2. Time-Integrated Extrapolated Fit 

An improved extrapolated fitting technique was devised such 
that a Band function was fit in the energy range from the vari- 
able LET to ~ 40 MeV. The function obtained from applying the 
model in this narrower range was then extrapolated to the lower 
energies in order to ascertain how well the function described the 
data below the LET. Deviations present between ~ 8 ke V and the 
LET are quantified by summing the residuals between ~ 8 keV 
and the LET. 

One of the assumptions of the extrapolated fitting technique 
is that the function parameters obtained over a shorter energy 



range should be consistent with a fit performed over the entire 
energy range in the absence of any spectral excesses or deficits 
with respect to the function. The parameters obtained by fitting 
from 8 keV - 40 MeV should be consistent with those obtained 
by fitting from 15 keV - 40 MeV when no additional compo- 
nents are present. The difference between the single fit method 
and extrapolated fit method is shown in Figure [3] Minimizing C- 
stat will tend to produce a fit that balances excesses and deficits 
regardless of the physical origin of the deviation from a Band 
function. In the single fit method, the deviation at low energies 
results in strong fluctuations throughout the entire spectral en- 
ergy range. When the extrapolated fit method is used, the spec- 
tral deviations are observed more clearly at lower energies and 
the effects on the remaining spectrum are reduced. Features in 
certain energy bands are not physically comparable between dif- 
ferent GRBs because the energy bands are defined in the ob- 
server frame and redshifts are not known for all GRBs. However, 
when performing an analysis on the manifestation of observable 
deviations in certain energy bands, fits using the same LET can 
be compared using this method. 

The low -energy power-law index, a, must be reasonably well 
constrained in order to confidently calculate the significance of 
any deviations. This implies that £ pea k an d the LET must be suf- 
ficiently different to ensure that a can be constrained over the 
shortened energy range. The technique also requires that £p ea k is 
greater than the LET. 

For each GRB, the spectral residuals between 8 keV and the 
LET are summed in an individual detector. The overall value 
for a particular GRB is the arithmetic mean for deviations in all 
detectors used in the analysis of the burst. The Iodine K-edge 
region of the spectrum around 33 keV can lead to larger resid- 
uals caused by systematics rather than a process intrinsic to the 
source. However, a strong K-edge effect was only noted in ~ 
5 GRBs and dominated only 2 spectral bins. If the K-edge is 
present between 8 keV and the LET, it only weakly affects the 
summing of residuals over this range. In reality it is only rele- 
vant when LET = 50 or 100 keV and the effect was not deemed 
significant enough to alter the method. 
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Fig. 3. Comparing the single time-integrated fit to the extrapolated fitting technique, a) Single Band function fit to GRB 090902B 



from 8 keV - 40 MeV (a = -0.99!;;{j;, E 
MeV and extrapolated down to 8 keV (a 
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3.3. Time-Resolved Extrapolated Fit 

The spectrum of a GRB can evolve over the prompt emission 
interval (e.g. |Preece et al.|2000 1 and a time-resolved analysis is 
required to account for such evolution. An S/N approach was 
employed over a basic time-resolved analysis (i.e. splitting the 
bursts into time sections) to ensure significant counts per bin. 
The selected region for each GRB used in the time-integrated 
fitting (full burst) was binned to the 25<x and 50cr level above 
background in the brightest detector. These intervals were then 
used to define bins for all other detectors in a particular GRB. A 
fit was then performed on each newly defined bin. 



3.4. Simulations 

The goodness of fit of the data below the LET and the function 
fit at higher energies cannot easily be evaluated using C-Stat as 
there is no analytic result available to convert from the fit statistic 
to a measure of goodness of fit. Therefore it is difficult to quan- 
tify the significance of a deviation or even an acceptable range of 
values which occur in the absence of an excess/deficit at low en- 
ergies. Simulations are therefore necessary to quantify whether 
a deviation is consistent with the extrapolated fit. 

A boot-strapping method was used to compare the actual re- 
sults from the data with simulated results. In order to perform 
the simulations, the background and source must be simulated. 
The background was simulated using a similar level to that in 
the real data. To obtain the source region, a perfect Band func- 
tion specific to each GRB, or time interval selected, was simu- 
lated. Multiple count distributions were created using this func- 
tion (varying by Poisson noise). Each set of parameters was sim- 
ulated ~ 1000 times. 

For each simulated spectrum, a distribution of low-energy 
deviations using different LETs can be compiled. The simulated 
distributions were then fit with a Gaussian distribution to obtain 
a value for the mean and standard deviation (see Figure [5j. The 
data obtained by summing the low-energy residuals below the 



LET were then normalised by calculating the number of standard 
deviations by which the data varied from the mean of the simu- 
lated distribution. If the data for a GRB differ significantly from 
its simulated distribution, then it is claimed that a low-energy 
deviation is present. 

A representation of the overall distribution was simulated 
to provide an initial insight into the expected distribution of 
residuals assuming that all spectra were correctly described by 
a Band function. Five GRBs were selected from the sample 
which represented a broad range of peak energies, with £p ea k 
= 161, 276, 444, 484 and 1010 keV. These GRBs were simu- 
lated and the resulting distribution of residuals in the absence 
of excesses/deficits are presented in Figure [4] Individual simu- 
lations were then performed on a range of time-integrated and 
time-resolved sections of GRBs. 



3.5. Data Cuts 

Several cuts were applied to the data before analysis. If any fit 
failed in the automated fitting process due to a lack of spectral 
constraints, it was excluded from the sample. A range of initial 
parameters were passed to the fitting process in order to obtain a 
satisfactory solution, but 9 fits failed in the time-resolved fitting 
process despite alternative initial parameters and were excluded 
from the sample. These fits were individually examined and the 
failure was usually due to low £p ea k and a lack of counts above 
£peak- No fits failed in the time-integrated analysis. 

As the low-energy spectral index a only approaches a power- 
law in the asymptotic limit, a sufficient data range is needed to 
obtain reasonable constraints on the index. If E pes k is too close 
to the LET, two issues may arise 1) a will not approach the limit 
and 2) there will not be sufficient data to constrain a resulting in 
large fluctuations (e.g. LET = 50 keV, £p e ak = 60 keV). A well 
constrained a parameter is required because slight variations in 
a can cause large spectral deviations when extrapolated to lower 
energies. For a spectrum to be accepted in the sample the follow- 
ing criteria must be met: 
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Fig. 4. Comparing the time-integrated data to a simulated distribution, a) Combined distribution of low-energy deviations of 5 
individual GRBs assuming a perfect Band function (no deviations present). The extrapolated fit method was used with LET = 25 
keV. b) The low-energy residuals of the GRBs in the sample that survived the data cuts (see § 3.5 1 for LET = 25 keV. The data 
distribution b) is broadly similar to the simulated distribution a). The data point outside the main data distribution (on the right) is 
GRB090902B. 



Table 2. Number of sample GRBs / GRB intervals pre- and post- 
cuts, where TI are the time-integrated results and 25<x/50cr are 
the time-resolved results. 



LET 


TI 


TI 


25tr 


25o- 


50<x 


50cr 




Pre 


Post 


Pre 


Post 


Pre 


Post 


15 


45 


42 


672 


457 


309 


230 


20 


45 


41 


672 


438 


309 


227 


25 


45 


41 


672 


389 


309 


222 


30 


45 


41 


672 


338 


309 


217 


50 


45 


23 


672 


189 


309 


117 


100 


45 


22 


672 


81 


309 


83 



- For LET = 15, 20, 25, 30 keV, E peak > 100 keV is required 

- For LETs = 50 or 100 keV, E pea t, > 200 keV is required 

- All intervals with a large proportional error on E pea k, A 
£peak/£peak > 0.45 were excluded. 

- All intervals with an error on a > 0.2 were discarded to en- 
sure that a can be reasonably extrapolated beyond the range 
of the fit. 

The final sample after the cuts contains GRBs and intervals with 
good statistics and well constrained spectral parameters. Table 
|2] shows the number of GRB intervals both before and after the 
cuts were applied for the time-integrated and time-resolved anal- 
ysis. 

3.6. Variance of a 

Another indication that a time interval in a burst had a low- 
energy spectral deviation was a variation in a when only the 
LET is changed. For each time section, each a parameter ob- 
tained from using different LETs was compared. If the value of 
a remains consistent across all different LETs in a spectrum, it 
is assumed that no significant low-energy deviations are present. 
If variation in a of > 2 <x for different LETs was evident, the 



GRB was investigated further. Usually a burst with a strong de- 
viation in a would vary across multiple different LETs. If an 
interval had a variation in a and there was a deviation present 
in the low-energy residuals, this interval was concluded to have 
a significant spectral deviation from a Band function. All GRBs 
that had large deviations in the low-energy residuals displayed 
the expected variations in a. 

3. 7. Summary of Method 

A brief summary of the method for time-resolved analysis is pre- 
sented below. 

1. The highest fluence GRBs from the first 2 years were se- 
lected for analysis. 

2. For each burst, Nal detectors < 60° to the source and without 
blockages were selected. Multiple Nals and a single BGO 
were selected in most cases. 

3. Background and source regions were selected for an individ- 
ual burst. 

4. The main emission (time-integrated) interval of the GRB 
was selected by eye. 

5. The GRBs were binned by S/N over the time-integrated in- 
terval of the GRB to produce the time-resolved intervals. 

6. The individually significant S/N bins were then fit with a 
Band function from a Low-Energy Threshold (LET) to ~ 40 
MeV. 

7. The LETs were selected to be 15, 20, 25, 30, 50, 100 keV. 

8. The fit function was then extrapolated to 8 keV. 

9. The residuals were summed between 8 keV and the LET for 
each detector and averaged across detectors. 

10. Simulations were then run and fit with a Gaussian distribu- 
tion. 

1 1 . The significance of deviations in the data were then quanti- 
fied by the number of standard deviations by which the data 
varied from the mean of the simulated distribution. 
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12. The a parameters between different LETs for a time region 
were compared. 

13. If a time region had a significant residual deviation and a 
variation in a while only changing the LET, it was consid- 
ered to be a spectral deviation. 



4. Results 

In the context of this work, spectral deviations can be classified 
into two broad categories - excesses and deficits. The former 
occur when there is an excess of photons with respect to the ex- 
trapolated fit function below the LET. These, for example, could 
be explained by an additional component such as a power-law or 
blackbody dominating at low energies. Deficits occur, not due to 
a lack of photons at low energies but because a is flatter than the 
actual data at low energies. This can be caused by the presence 
of an additional component, such as a blackbody, between the 
LET and £p ea k resulting in a shallower a and a deficit of photons 
relative to the extrapolated function at low energies. 

4. 1 . Hypothesis Testing 

Two methods were used to test for deviations at low energies. 
The first involves comparing the residuals from the data to those 
from simulations where a Band function is the assumed best fit. 
The results in Figure[5]show two GRBs, one with and one with- 
out deviations using the time-integrated results. The null hypoth- 
esis for these simulations is that if the data are consistent with 
the simulated distribution, the extrapolated spectrum is consis- 
tent with a Band function. 

The second method involves looking at the variation in a 
while only changing the LET for the time-resolved analysis. The 
changes in a for one GRB which is consistent with a Band func- 
tion and one which is not consistent are presented in Figure [6] 
Two GRBs, GRB 090618 and GRB 091024 showed variations in 
a but did not have large deviations in the low-energy residuals. 
The parameter a varied outside 2 cr during one temporal interval 
for each of the GRBs. These intervals were investigated further 
and were discounted due to a long, weak interval in GRB 091024 
and only having a single 'good' detector available for analysis in 
GRB 090618. 

Simulations were performed for GRBs that initially pre- 
sented strong deviations in the time-integrated results. These de- 
viations were then quantified by comparing the data to the sim- 
ulated distribution for that GRB and were discarded as a chance 
occurrence if consistent with the simulations. Otherwise it was 
classified as a spectral deviation. A sample of GRBs were sim- 
ulated that did not show strong deviations to ensure the veracity 
of the null hypothesis. 

Simulations could not be performed for time-resolved re- 
gions of all GRBs so the combination of the two methods out- 
lined previously (summing the low-energy residuals and check- 
ing the change in a between different LETs) present the final 
criterion by which a deviation is tested. If a deviation passed 
both tests then it was classified as a real spectral anomaly. 

4.2. Time-Integrated Spectra 

Out of the sample of 45 GRBs tested, a significant devi- 
ation in the time-integrated results was found in 3 cases 
GRB 090902B (T0+0 - TO+25 s), GRB 090424 (T0+0 - TO+6.4 



s), GRB081215A (T0+0 - TO+8.4 s). A large 28.2 cr excess 
(normalised) is noted in GRB090902B and is a clear outlier 
in Figure |4] GRB 090424 has a deficit of 6.1 cr in the time- 
integrated spectrum which appears to be due to an anomalous 
spectral feature near Epeak- These two GRBs also show devia- 
tions in the time-resolved spectral analysis. GRB 081215A man- 
ifests as a 6.7 cr deficit but does not show up significantly in 
the time-resolved analysis. This deficit appears to be caused by 
strong spectral evolution throughout the GRB with £p ea k de- 
creasing from > 1 MeV in the first 1 s time bin to ~ 50 keV 
at the end of the 8 s emission region (a time-resolved analysis 
by |Zhang et al. ( 2011| l using different binning shows similar re- 
sults). 

GRB 090902B may be fit with an intense additional power- 
law component throughout the entire burst causing it to devi- 
ate strongly in the time-integrated spectrum ( |Abdo et al.|2 009a 



Zha ng et al.|2011[ ). GRB 090424 has an intense additional com 
ponent, which can be modelled as a blackbody throughout the 
burst which also manifests as a deviation in the time-integrated 
results. Additional emission processes and components may al- 
ways be present in all GRBs but the instrument may not be sen- 
sitive enough to confidently detect them. A closer inspection of 
each GRB using time-resolved analysis is required to differen- 
tiate between artifacts in the spectra and real additional time- 
dependent spectral features. 



4.3. Time-Resolved Spectra 

A time-resolved analysis was performed on 672 spectra with an 
S/N of 25 cr and 309 spectra with an S/N of 50 cr which resulted 
in significant features in five bursts. GRB 090323, GRB 090424, 
GRB 090820, GRB 090902B and GRB 090926A contained at 
least one temporal interval where significant deviations were de- 
tected. Simulations were performed on the relevant time regions 
to ensure that they were significant. The properties of the five 
bursts are presented in Table [3] 

Out of the five bursts observed with low-energy deviations, 
four are in the top five most fluent GRBs from the first 2 years. In 
the case of the highest fluence event ( GRB 090618 with a fl uence 
of 2.68 x 10~ 4 +4.29x 10~ 7 erg / cm 2 ( jPaciesas etal.|2012| », only 
one Nal was available for analysis due to the source-instrument 
geometry. It cannot be claimed with certainty that the result is 
stable on the basis of one 'good' Nal detector. An analysis of this 
burst with the combined GBM-Swift data did not significantly 
demonstrate the need for an additional component in the prompt 
gamma-ray emission (Pa ge et al.|2011| l. 

As spectral evolution can lead to anomalous behaviour, even 
in the time-resolved analysis, lightcurves were examined for any 
notable features during a time interval which contained a devia- 
tion. To limit the effect of temporal binning, several further fits 
were performed by changing the length of the interval by up to + 
50 % and by shifting the interval up to 50 %. If the effects were 
still significant after this further analysis it is claimed as a detec- 
tion of an additional spectral feature. These intervals were then 
fit with more complex models and the results are summarised 
in Table |4] There is a notable improvement in the fit statistic 
when a Band function is fit with an additional component. This 
provides evidence that additional spectral features are needed to 
accurately model these time intervals. These solutions however 
are not unique and other spectral models may have a similar or 
better fit statistic. 
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Fig. 5. Histograms of summed residuals below 25 keV obtained from time-integrated simulations of a perfect Band function. The 
dark line shows the value of the data for that GRB. a) The value of the deviation for GRB 0808 17A is -3.6. The mean and the 
standard deviation of the simulated distribution are -1 . 1 and 3.4 respectively. In this case the data does not deviate significantly from 
the simulated distribution and has a normalised deviation of -0.7 cr. b) The value of the deviation for GRB 090424 is -25.8 compared 
to the simulated distribution with a mean and standard deviation of 0.4 and 4.3 respectively. The data deviates significantly from the 
simulated distribution with a normalised deviation of -6.1 cr. 

Table 3. Properties of the GRBs found to have strong spectral deviations in the time-resolved analysis at low energies. 



GRB 


Fluence" 




Fluence Rank" 


T90 fl 


Deviation'' 


Localisation 1 ^ 




10 - 1000 keV 


first 2 years 


50 - 300 keV 








(erg / cm 2 ) 






(s) 






GRB 090323 


1.2 x 10" 4 


± 1.7 x 10- 7 


5 


135.2+1.5 


+4.9 cr (Excess) 


Swift XRT 1 


GRB 090424 


4.6 x 10~ 5 


± 3.9 x 10- 8 


19 


14.1±0.3 


-6.2 cr (Deficit) 


Swift XRT 2 


GRB 090820 


1.5 x 10~ 4 


± 1.8 x 10~ 7 


3 


12.4±0.2 


-5.5 cr (Deficit) 


GBM 3 


GRB 090902B 


2.2 x 10~ 4 


± 3.2 x 10~ 7 


2 


19.3±0.3 


+25.0 cr (Excess) 


Swift XRT 4 


GRB 090926A 


1.5 x 10~ 4 


±3.4x 10~ 7 


4 


13.8±0.3 


+6.5 cr (Excess) 


Swift XRT 5 



Notes, (a) From Paciesas et al. (2012). (b) The deviations from the extrapolated fitting method have been normalised by the standard deviation 
from the simulated distributions for each individual burst section, (c) The Swift XRT location is sufficiently accurate to model the response of the 
instrument for spectral analysis, even in cases where a more accurate location is known. 

References. (l)|Kennea eUtf"| (20091); (2)|CannizzolFaT]((2009)l; (3)|Connaughton (2009b); (4)|Kennea & Stratta (2009 ); (5)|Vetere et aT"](|2009) 



4.4. Individual Bursts with anomalous low-energy spectra 

Significant anomalous spectral behaviour with respect to a Band 
function was found in time intervals of five bursts. The spectra 
were further investigated by fitting a number of spectral models 
as presented in Table|4] The individual events are described in 
more detail below. 



4.4.1. GRB090323 



2009 ). The energetics and host galaxy properties of this burst are 
described in |McBreen et al.lpOTO) ). 

After performing a spectral analysis on this burst an excess is 
detected in two out of nine intervals (52 - 60 s, 60 - 66 s relative 
the trigger time). The spectral analysis for the interval with the 
strongest low -energy deviation is presented in Figure|7]and Table 
[4] An improvement in the spectral residuals is notable when a 
Band + blackbody or Band + power-law fit is performed over a 
single Band fit. 



GR B 0903 23 was detected by Fermi GBM ( |van der HorsT&l 
Xin 2009) and LAT which caused the satellite to perform an 
autonomous repoint to the source location (Ohno et al. 2009). 



The lightcurve from GBM Nal detector n9 is displayed in Fig. [7] 
and consists of multiple distinct pulses over ~ 150 s. Multi- 
wavelength observations were aided by an on-ground location 
provided by the LAT and X-ray observations from Swift-XRT 
( Kennea et al.||200"9] ). Follow-up observations report a spectro- 
scopic redshift of z = 3.57 ( |Chornock~e t al. 200 9a]) and the af- 
terglow was also detected at radio wavelengths ( [Harrison et ak] 



4.4.2. GRB 090424 

The prompt emission from GRB 090424 was detected by Swift- 
BAT, XRT and UVOT which provided an arcsecond localisa- 
tion to the follow-up community within minutes (Cannizzo et al. 
2009). The prompt emission was also observed by Fermi GBM 
( |Connaughton|2009a| ) and Suzaku WAM ( |Hanabata et al.|2009 1. 
Multiple optical telescopes observed the optical afterglow (e.g. 
( jOlivares et al.|2009b| )) and a spectroscopic redshift of z = 0.544 
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Fig. 6. Change in a between different LETs. The S/N temporal selections in GRB 090102 and GRB 090902B are presented in a) 
and c). The spectral parameter a is presented with 2 <x error bars in b) and d) for LET values from 15 keV to 50 keV. a is consistent 
at the 2 cr level as the LET is varied for the time interval selected in GRB 090102 and is not consistent for the time interval in 
GRB 090902B. 



was obtained by Gemini-South (Chornock et al. 2009b). Radio 



observations detected a bright radio afterglow ( Chandra & Frail 
2009) 1 and a disturbance in Earth's ionosphere was detected by 



monitoring very low frequency radio waves ( |Mondal et al. 201 1 
IChakrabarti et al.|2010) . 

The main emission period of this GRB is relatively short 
compared to other bursts in the high fiuence sample and has 



an extremely high peak flux of 110 photons / s (Paciesas et al. 
|2012| . The burst has at least 3 pulses in rapid succession last- 
ing ~ 6 s and lower level emission that continues up to 20 s 
after the trigger time. A significant low-energy deficit was ob- 
served in the time-integrated and time-resolved spectral analysis 
of GRB 090424. This apparent deficit is caused by the unusual 
spectral shape, which is shown in Figure [8] The spectrum has 
2 distinct spectral breaks, one in the usual energy range where 
a break due to £p e ak is to be expected but also another lower 
energy break in the spectrum. A Band + blackbody fit with a 
blackbody kT of ~ 9 keV can be fit throughout the burst which 



is most prominent between 2 - 3 s relative to the trigger time. Six 
1 s intervals between - 6 s were analysed and a strong deviation 
is present in the residuals in the the 2 - 3 s region. Little spectral 
evolution is present throughout the burst, and as this feature is 
observable throughout the burst it is less likely to be caused by 
spectral evolution. 



4.4.3. GRB 090820 

GBM detected GRB 090820 and issued an automated repoint re- 
quest to Fermi, however the LAT could not observe the burst due 
to Earth avoidance constraints ( |Connaughton|200 9b ). The RT-2 
Experiment onboard the CO RONAS-PHOTON satel lite also ob- 
served the prompt emission (Chakrabarti et al.|2009) . No further 
follow-up of the burst occurred and thus the best location for 
the burst is obtained from GBM. The burst triggered on a weak 
precursor which was excluded from the spectral analysis. 
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Table 4. Time-resolved spectral parameters from model fits to GRB time intervals where significant deviations were detected. 
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Notes. *) Measured in seconds since trigger time (TO). 1) The Band+PL fit for GRB 090820 did not converge to physical parameters and is 
omitted. 2) /3 was frozen to -3.00 if the uncertainties were unconstrained. 3) Although many intervals have significant excesses in GRB 090902B, 
only one is presented above for illustrative purposes (See also Abdo et al. 2009a, P e'er et aTf2 012; Ryd e et al.|2011) . 



This burst presents a low-energy deficit relative to a Band 
function which is caused by the Band function attempting to fit 
data that has a broader peak relative to the expected function fit. 
Thirteen time intervals between 29 - 45 s from the trigger time 
were analysed with five 1 s intervals between 31 - 36 s exhibit- 
ing marginal evidence of deviations. As this was the only GRB 
with consecutive marginal evidence for deviations, these time in- 
tervals were analysed together and display a strong deviation in 
the spectral data. By adding a blackbody component the Band 
function can account for the broader spectral peak as shown in 
Figure [9] As an additional power-law could not account for de- 
viations at lower energies, this provides evidence that some ad- 
ditional spectral property is occurring in the mid-energy range. 
Spectral analysis comparing standard models to more advanced 
fitting techniques, such as synchrotron models have been carried 
out by Burgess et al.| ( |20lT[ ) for this burst. 



4.4.4. GRB090902B 



GRB090902B (Bissaldi & Connaughton 2009) was an ex- 
tremely bright GRB with a bright additional component ob- 
served across 8 orders of magnitude in energy (Abdo et al. 
2009a). The GRB was initially detected by both instruments 
on Fermi, GBM and the LAT ( |Bissaldi & Connaughtonl|2009| 
de Palma et al. 200 9~]l. The prompt em ission was also detected 



by Suzaku WAM ( |Terada et al.|200 9). A Target of Opportunity 
(ToO) was issued to Swift which observed an uncataloged source 
within the LAT error radius (Kennea & Stratta 2009). Numerous 



optical follow-up observations were made (e. g. |Qlivares et aL 
2009a) resulting in a redshift of z = 1.822 ( |Cucchiara et al. 
2009). The source was also observed in the radio ( |van der Horst 



et al. 2009). 



A significant excess was observed in this GRB in numerous 
intervals (ten 1 -second intervals between 6 - 16 s relative the 
trigger time). The region with the stongest deviation is presented 



in Table [4] An additional power-law component was the best fit 
model to the data out of the models tested. The additional power- 
law index from analysing the GBM data only in this region, a = 
-2.00!jjgg, is consistent with the power-law index observed by a 
joint fit over the entire spectral range of GBM and LAT {a - - 
'-0 02^ Complex spectral models to this GRB are fit in Figure 
Numerous models have been presented in the literature to 



10 



explain the spectral features in this burst (Liu & Wang 2011 



B ar niol Duran & Kumar|20 1 1 [ IPe^er et al. |20 1 2) . 



4.4.5. GRB090926A 



Detec ted by GBM ( |Bissaldi| [2009) i and LAT ( |Uehara et aL] 
|2009[ ), follow-up observations of this GRB were made by Swift- 
XRT HVetere et al.|2009) and Swift-UVOT ( |Gronwall & Vetere| 



2009). The prompt emission was also observed by Suzaku WAM 
flNoda et al.l|2009l and Konus-Wind ( |Golenetskii et al.||2009) . 
A redshift of z = 2.1062 was obtained by the Very Large 
Telescope (VLT) X-Shooter spectrograph (Malesani et al. 2009|l. 
Skynet/PROMPT also observed the optical afterglow (|Haislip 
|et al.|2009D . 

GRB 090926A has also been observed in joint GBM and 
LAT spectral fitting to necessitate an additional power-law com 
ponent over a wide spectral range (Ackermann et al. 2011 



Thirteen time bins were analysed between - 17 s with one time 
bin showing a significant excess. The excess is observed during 



a short sharp spike in the lightcurve which is shown in Figure 1 1 



When an additional power-law is fit with a Band function, the fit 
statistic shows a large improvement. The additional power-law 



index a - -1.73*QgJ, obtained by using GBM data only, is again 
ower-law observed in t he LAT at high ener- 



-0.04 

consistent with the j: 
giesa = -1.71!°°? 



Ackermann et al 



2011 
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Fig. 7. Spectral fits for GRB 090323. a) A Band spectral fit and b) a Band + blackbody spectral fit are displayed for the hatched 
region of interest in c) the CSPEC lightcurve. All fits were performed on the full resolution data and the data were binned for 
illustrative purposes. The fit parameters are given in Table |4] 



5. Discussion and Conclusions 

Although an empirical model, a Band function has been used to 
model the majority of GRB spectra since it was first proposed 
in 1993 ( [Band et af|l993) . However recent GRB studies with 
Fermi have suggested that there are a number of GRBs whose 
spectral behaviour differs from a Band function, with a num- 
ber of different spectral components invoked in order to explain 
these observations (e.g. Asano et al.|20 09 ; Burg ess et al.|20ir 



Guiri ec et al.||20lf] |Pe'er et al.||2012| l. The Fermi data cover a 
wide energy band with high spectral resolution and so the sys- 
tematic study of low-energy spectral deviations in bright Fermi 
bursts presented here enables the frequency of such events to be 
investigated. Spectral deviations to a Band function occurring 
throughout the GBM energy range also affect the spectrum at 
low energies, and so the current method provides a systematic 
study over the full GBM energy range, assuming that the spec- 
tral data are adequately modelled by a Band function. 

GBM observes 3 out of 45 (~ 7%) GRBs with spectral de- 
viations present in the time-integrated data. The probability that 
our observed rate is consistent with the BATSE observed rate 



is ~ 7.0% or within 1.8 <x. An excess is observed in one GRB 
(GRB 090902B) and deficits are observed in two (GRB 090424 
and GRB081215A) in the time-integrated data, while BATSE 
observed only excesses at low energies. There are a number of 
possible causes for the differences between the results obtained 
from the two instruments. As described in § 1.2 the study of 
low-energy GRB spectra with BATSE was performed using two 
data types from two different SDs, with the observations de- 
pending on a single DISCSP1 data point at low energies. The 
structure of these datatypes mean that the low and high-energy 
bounds of the DISCSP1 point can vary from one GRB to an- 
other. GBM has the advantage of a single data type with high 
spectral resolution over its full energy range allowing deviations 
to be investigated in fixed energy ranges. Multiple GBM detec- 
tors were also used for this analysis, with GRBs having only one 
'good' detector being analysed for 3 out of 45 GRBs. The effec- 
tive areas of both instruments are also different. The DRMs of 
the GBM GRBs are well defined due to their small source loca- 
tion uncertainties, with 80 % of the GRBs with spectral devia- 
tions having Swift-XRT locations. It is therefore difficult to make 
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Fig. 8. Spectal fits for GRB 090424. a) The high-time resolution lightcurve of GRB 090424 with the hatched region being the region 
of interest, b) A simple Band fit to selected region, c) A Band + blackbody fit to the data, d) A double broken power-law fit with 
parameters of E breakl = 32.74! J-g keV, E break2 = 197.70^ |° keV, Indexl = -0.59^, Index2 = -1.52+^, Index3 = -2.95!°;}* 
(where Indexl < E break i, E blea ki < Index2 < E break2 , Index3 > E break2 ). 



a direct comparison between the GBM and BATSE results. For 
the reasons outlined above we consider that the GBM result is 
more robust. 

Since spectral evolution throughout the duration of a GRB 
significantly affects the spectral shape, a time-resolved analysis 
was performed for all of the GRBs in the sample. Spectral devi- 
ations are found in 1 1 % of the GRBs in at least one time interval 
by applying the methods outlined in this paper. Our work shows 
that the features in the spectra are time-dependent and a time- 
resolved analysis is necessary to limit the effects of averaging 
multiple spectra. It is interesting to note that 4 of the 5 GRBs 
with significant features in the time-resolved analysis are in the 
top 5 most fluent GRBs in the sample. A bias exists whereby 
additional models are more likely to be required in high fluence 
GRBs since these GRBs have enough statistics to confidently 
define an additional component, if present. 

The observed deviations may be due to significant spec- 
tral evolution in the burst (or time interval) or additional spec- 
tral features/components in the spectrum. In the time-integrated 



analysis, the deviation found in GRB 08 12 15 A may be due to 
spectral evolution, while the deviations in GRB 090424 and 
GRB 090902B could be explained by additional spectral fea- 
tures. For those GRBs with deviations found in shorter time 
intervals (GRB 090902B and GRB 090926A), a power-law im- 
proves the fit at lower energies which is consistent with the 
high-energy observations in the LAT data (Ab do et al.| [2009a; 
|Ackermannetal.|2011) . For GRB 090424 and GRB 090820 an 
additional blackbody component provides a significantly better 
fit to the data than an additional power-law, showing that a num- 
ber of different spectral components may be present in GRBs. 
GRB 090323 is adequately fit invoking either of the two addi- 
tional components. However as there is little evidence for LAT 
emission ( Piron et al.||201 l| l, it is unlikely to be due to an ad- 
ditional power-law component extending to high energies and 
a spectral cut-off between the GBM and LAT energy bands is 
possible. 

The advantage of the method presented in this paper is that it 
can demonstrate the requirement for an additional spectral com- 



12 



D. Tierney et al.: Anomalies in low-energy GRB spectra with Fermi GBM 




ponent with respect to a Band function without any prior knowl- 
edge of the nature of that extra component. Theoretical models 
of the GRB emission process must be able to explain why the 
majority of GRBs are adequately modelled by a Band function 
and also explain the observed deviations. Many physical models 
have recently been proposed (see § [TJ, however in many cases 
the spectral data do not have adequate statistics for one model to 
be deemed more statistically significant than another. 

This technique demonstrates a systematic method to search 
GRBs for additional components at low-energies. The method 
finds several GRBs in which features have previously been 
noted as well as discovering deviations in GRB 090424 and 
GRB 090323. Deviations are found in a large fraction of high 
fiuence GRBs; this method may be unable to find low-energy 
deviations in fainter GRBs. This suggests that deviations from a 
Band function may be common in all GRBs but the statistics do 
not allow us to test this hypothesis. 
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